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ABSTRACT 

Using deep Keck spectroscopy of Lyman break galaxies selected from infrared imaging data taken 
with the Wide Field Camera 3 onboard the Hubble Space Telescope, we present new evidence for a 
reversal in the redshift-dependent fraction of star forming galaxies with detectable Lyman alpha (Lya) 
emission in the redshift range 6.3 < z < 8.8. Our earlier surveys with the DEIMOS spectrograph 
demonstrated a significant increase with redshift in the fraction of line emitting galaxies over the 
interval 4 < z < 6, particularly for intrinsically faint systems which dominate the luminosity density. 
Using the longer wavelength sensitivities of LRIS and NIRSPEC, we have targeted 19 Lyman break 
galaxies selected using recent WFC3/IR data whose photometric redshifts are in the range 6.3 < z < 
8.8 and which span a wide range of intrinsic luminosities. Our spectroscopic exposures typically reach 
a 5cr sensitivity of < 50 A for the rest-frame equivalent width (EW) of Lya emission. Despite the 
high fraction of emitters seen only a few hundred million years later, we find only 2 convincing and 
1 possible line emitter in our more distant sample. Combining with published data on a further 7 
sources obtained using F0RS2 on the ESO VLT, and assuming continuity in the trends found at 
lower redshift, we discuss the significance of this apparent reversal in the redshift-dependent Lya 
fraction in the context of our range in continuum luminosity. Assuming all the targeted sources are 
at their photometric redshift and our assumptions about the Lya EW distribution are correct, we 
would expect to find so few emitters in less than 1% of the realizations drawn from our lower redshift 
samples. Our new results provide further support for the suggestion that, at the redshifts now being 
probed spectroscopically, we are entering the era where the intergalactic medium is partially neutral. 
With the arrival of more sensitive multi-slit infrared spectrographs, the prospects for improving the 
statistical validity of this result are promising. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: starburst - galaxies: high 
redshift 



L INTRODUCTION 

Determining when neutral hydrogen in the intergalac- 
tic medium (IGM) was reionized is an important ques- 
tion in observational cosmology and a precursor to un- 
derstanding whether star forming galaxies prov ided the 
necessary UV photons (jRobertson et al.l l2010f ). One 
of the most practical probes of reionization with cur- 
rent facilities utilizes the frequency of occurrence of 
Lya emission in star forming galaxies. As Lya pho- 
tons are resonantly scattered by neutral hydrogen, the 
abundance of Lya emitters should decrease as ob- 
servations probe into the era whe r e there is neutral 
gas (e.g., 'Malhotra & Rhoadd [IQOI lOuchi et al.l [2010l 
[Kashikawa et al. 2011). The recent discovery of large 
numbers of candidate galaxies beyond z ~ 7 through 
multi-color imaging undertaken with the infrared Wide 
Field Camera (WF C3/IR) onboard Hubb l e Space Tele- 
scope (HST) (e.g., iBouwens et"all l2010bt iBunker et all 
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[20T0I : iFinkelstein eTaH IMot iMcLure et all [20To[ ) now 
makes it feasible to track the occurrence of Lya line emis- 
sion to interesting redshifts where neutral hydrogen may 
be present. 

Of course, astrophysical factors other than a neutral 
IGM can also affect the presence of Lya emission. Be- 
cause of this, an alternative approac h for gaugin g when 
reionization occurred, introduced by iStark et all ()2010f ) 
(hereafter Paper I), is to spectroscopically measure the 
fraction of Lya emission within color- selected Lyman 
Break Galaxy (LBG) populations. Tracking the redshift- 
dependent fraction in a well-defined population avoids 
consideration of absolute changes in the host galaxy num- 
ber density such as has been the case in studies based on 
the luminosity function of narrow-band selected Lyman a 
emitters (LAEs). Furthermore, evolution in dust obscu- 
ration can be independently tracked via correlations seen 
with the colors of the rest-frame UV continuum. Paper I 
presented a comprehensive survey of over 600 LBGs with 
deep spectra, mostly undertaken with the DEIMOS in- 
strument on Keck, but including pub lished sampl e s from 
the Very Large Telescope (VLT, iVanzella et al.l (|2009f ) 
and references therein). In that paper we demonstrated 
the utility of the method and discussed the paucity of 
line emission in gr avitationally-lensed z > 7 candidates 
from the sample o f iRichard et al.l (|2008D . 

In IStark et al.l (J2011D (hereafter Paper II), through 
ultra-deep exposures with DEIMOS we significantly im- 
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proved the line emission statistics at z ~ 6, providing a 
robust measure of the rest-frame EW distribution of Lya 
emission at the highest re dshift when the U niv erse is be- 
lieved to be fully ionized (jFan et alJ 120061 c.f. iMesingeJ 
I201CO . This provides a sound basis for predicting the like- 
lihood of emission at higher redshift and thereby enabling 
a test of whether there is absorption by neutral gas. Sig- 
nificantly, we found that over 50% of moderately-faint 
(-20.25 < Muv < -18.75) z ~ 6 LBGs exhibit strong 
emission with rest frame EWs > 25 A. As this fraction in- 
creases over 4 < z < 6, we argued on continuity grounds 
that we should expect a high success rate in recovering 
line emission from the newly-found WFC3/IR samples of 
z > 7 LBGs unless we encounter a more neutral IGM in 
the short time interval prior to z ~ 7. Some evidence for 
this is seen in the r e cent s tudies of Fontana et al. (2010) 
and lVanzeUaet al.l (|2011D . 

The present paper is concerned with an initial ap- 
plication of this test to the newly-available sample of 
WFC3/IR candidates with photometric redshifts in the 
redshift range 6.3 < z < 8.8. Two important factors 
have motivated and shaped our program. Firstly, it is 
important to note that Lya emission is the only spectro- 
scopic redshift indicator for galaxies beyond z ~ 6. Since 
it is the absence of strong Lya emission that provides the 
basis for considering an increased neutral fraction, it is 
important to be sure that the targets are truly at the 
expected redshifts. Many early candidate LBGs believed 
to lie beyond z ~ 6 — 7 remained controversial because of 
their limited or marginal photometry. The improved fil- 
ter set and superior performance of WFC3/IR has given 
us confidence that the current list of z ~ 7 — 8 candi- 
dates is mo re robust than those b ased on earlier NIC- 
MOS data (jRobertson et al.ll201Cil ). Secondly, to match 
the lower redshift data, sampling from a similarly wide 
range of LBG luminosities, as we do here, will be advan- 
tageous. As shown in Papers I and II, the fraction of line 
emission increases in intrinsically fainter sources and so 
by comparing fractions with respect to their LBG lumi- 
nosities, we may gain additional evidence for the onset 
of the neutral era. 

Throughout the paper, we adopt a A-dominated, 
flat universe with f^A = 0.7, Hm — 0.3 and Hq = 
70h7o kms^ Mpc~ . All magn itudes in this paper are 
quoted in the AI3 system (jOke fc Gunn 1983). 

2. OBSERVATIONS 

In compiling a target list for this program, we are 
guided by the need for a robust photometric redshift 
for each galaxy based on improved photometry from 
WFC3/IR and a range of rest-frame UV luminosities 
(Muvu)- Our primary source of targets for the wide- 
field multi-slit capabilities of the Low Resoluti on Imaging 
Spectrometer (LRIS) on the Keck I telescope (|Oke et al.l 
[1995:1 equipped with a new red-sensitive CCD was i' 
and z'-drop candidates whose photometric redshifts z > 
6.3 from the HST Early Release Sci ence (ERS) field 
(jHathi et al.ll2010HMcLure et al.ll2011h . The grating for 
these observations was blazed at 600 lines mm~^. On 
January 7 and February 4 2011 we secured 7 hours of 
on-source integration for 8 suitable targets on a single 

^ corresponding to a rest wavelength A ~ 1500 A (Paper I) 



mask using slit widths of 1", observed through a median 
seeing of 0.98". 

In a more ambitious campaign probing to higher red- 
shift we also targeted 3 z'-drop sources from the Hubble 
Uhradeep Field (HUDF) P34 field (GO 11563, PI: Ilhng- 
worth) and an additional gravitationally-lensed source 
in the cluster MS0451-03 (GO 11591, PI : Kneib) using 
the n ear-infrared spectrograph NIRSPEC (jMcLean et al.l 
[1991 during November 14-17 2010 and Jan 14-15 2011. 
This extends our search for Lya emission up to a red- 
shift z ~ 8.2. Although we undertook extended integra- 
tions on all 4 sources with a 0.76" slit, tracking difficul- 
ties affected some exposures. To determine the effective 
on-source integration time, we secured our astrometric 
position for each exposure by locating objects visible in 
the slit viewing camera to a precision of ^ 0.2". We 
continued this campaign over May 15-18 2011. During 
these 4 nights, we did not encounter any tracking dif- 
ficulties and, in excellent conditions, successfully used 
NIRSPEC to study an additional 7 WFC3-IR dropouts 
drawn fro m numerous survey s: the BoRG pure paral- 
lel survey (Trenti et al." '2011|) (also independently dis- 
covered bi^_Yan_rtal^01 1), the EGS region in CAN- 
DELS (|Grogin et all 120111 : iKoekemoer et all 1201 It ), and 
the lensing clusters Abell 1703 (GO 10325, PI: Ford, 
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20Tll) and Abeh 2261 (CLASH survey. 



We reduce the LRIS data following standard proce- 
dures, with bias subtraction and flat-fielding using dome 
exposures. We used the iKelson (2003:' ) code to remove 
spatial and spectral distortion and to model and sub- 
tract the sky emission. Wavelength calibration was de- 
termined directly from sky lines. A final two-dimensional 
spectrum was extracted for each object with pixels 
binned logarithmically by Alog(A) = 4.02 x 10~^. As in 
Papers I and II, we search through the two-dimensional 
spectrum visually to identify emission lines, and confirm 
these with a boxcar extracted one-dimensional spectrum. 

Our exposures with NIRSPEC were conducted with 
typical spatial dithering of 5". In the case of some lensed 
sources, we dithered by longer amounts to ensure the arc 
was oriented along the slit, and was not spatially coin- 
cident with any background objects. We flat-fielded and 
sky-subtracted the spectra using IDL routines written 
by G. Becker (2010, private communication). To com- 
pute the camera distortion and spectral curvature, we 
fit traces of standard stars along the slit, and skylines 
perpendicular to the slit. From this, we then derive a 
wavelength and sky position for all pixels in each two- 
dimensional spectrum, which were used to align the in- 
dividual exposures. A final spectrum was created for 
each object by median stacking all exposures to elimi- 
nate signals from cosmic ray strikes. 

In total, this paper therefore presents the results of 
Keck spectroscopy for 19 WFC3-IR selected sources 
whose photometric redshifts lie in the range 6.3 < z < 
8.8. A summary of the new observations is given in 
Table 1. To this sa mple, we add a furthe r 7 z > 6.3 
sources discussed by iFontana et al.l ()2010f ). Figure 1 
compares the UV absolute magnitude distribution of the 
combined sample with that presented for the redshift 
range z ~ 5 — 6 in Paper II; clearly the samples span 
a similar luminosity range. This luminosity range is 
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TABLE 1 

Catalog of sources with Keck spectroscopy and emission line properties, when detected, la magnitude errors are listed 

IN parenthesis. For non- detections, we list 2(t limiting magnitudes. For our NIRSPEC targets, times correspond to 

exposures in filters N1 and N2, which cover wavelength ranges of 9470 - 11210 A AND 10890 - 12930 A, RESPECTIVELY. 



ID 



R.A. 



Dec 



2850 



Jl25 



Hv 



^phot 



fhrl 



EW [A] 



LRIS 



ERS 5847^ '^ 


03:32:16.0 


-27:43:01.4 


26.6(0.1)^ 


26.6(0.1) 


26.7(0.1) 


- 


6.48 


7 


ERS 7376^ 


03:32:29.5 


-27:42:04.5 


27.2(0.1)1 


27.0(0.1) 


27.0(0.1) 


- 


6.79 


7 


ERS 7412 


03:32:10.0 


-27:43:24.0 


26.9(0.1)1 


27.0(0.1) 


26.7(0.1) 


- 


6.38 


7 


ERS 8119 


03:32:29.5 


-27:41:32.7 


27.7(0.2)1 


27.1(0.1) 


27.5(0.1) 


- 


6.78 


7 


ERS 8290 


03:32:13.4 


-27:42:30.9 


27.3(0.1)1 


27.1(0.1) 


26.8(0.1) 


- 


6.52 


7 


ERS 8496^* 


03:32:29.7 


-27:40:49.9 


27.2(0.1)1 


27.3(0.1) 


27.5(0.1) 


- 


6.52 


7 6.441 69 ± 10 


ERS 10270 


03:32:29.5 


-27:42:54.0 


28.1(0.1)1 


27.4(0.1) 


28.0(0.2) 


- 


7.02 


7 


ERS 10373 


03:32:27.0 


-27:41:42.9 


27.5(0.1)1 


27.4(0.1) 


27.8(0.2) 


- 


6.44 


7 


NIRSPEC 




A1703.zDl'= 


13:14:59.4 


51:50:00.8 


25.8(0.2) 


24.1(0.1) 


24.0(0.1) 


9.0 


6.75 


2. - 


A1703.zD3^ 


13:15:06.5 


51:49:18.0 


26.8(0.5) 


25.5(0.1) 


25.1(0.2) 


7.3 


6.89 


2,- 


A1703_zD6'= 


13:15:01.0 


51:50:04.3 


27.9(0.5) 


25.8(0.1) 


25.9(0.1) 


5.2 


7.02 


5, 3 7.045 65 ± 12 


A1703_zD7'' 


13:15:01.3 


51:50:06.1 


>28.5 


26.8(0.2) 


26.4(0.2) 


5.0 


8.80 


5, 3 


A2261.1 


17:22:28.7 


32:08:30.9 


>28.6 


26.9(0.1) 


27.3(0.1) 


3.5 


7.81 


5.7,- 


BoRG_58.1787.1420' 


14:36:50.6 


50:43:33.6 


> 27.91 


25.8(0.1) 


25.9(0.2) 


- 


8.27 


2, 3 


EGS.K1 


14:19:24.2 


52:46:36.2 


> 27.81 


25.3(0.1) 


25.4(0.1) 


- 


8.27 


2.5,- 


HUDF09.799*5 


03:33:09.1 


-27:51:55.4 


>29.1 


27.7(0.1) 


27.6(0.2) 


- 


6.88 


4.5,- 


HUDF09.1584«5 


03:33:03.8 


-27:51:20.4 


27.2(0.1) 


26.7(0.1) 


26.6(0.1) 


- 


7.17 


5.5,- 


HUDF09.1596 


03:33:03.8 


-27:51:19.6 


27.3(0.1) 


26.8(0.1) 


26.8(0.1) 


- 


7.45 


5.5,- 6.905? 30 ±15 


MS0451-03.10 


04:54:08.8 


-3:00:29.1 


>28.3 


26.7(0.1) 


26.9(0.1) 


50 


7.50 


2.5,- 



^Best fit magnification for sources in our sample wiiich are gravitationally lensed. 

''Y098A/ 

'^l814VK 

"^In order listed in table, discovered as Objects 1,4,3 in lWilkins et al] II2010I) 

"^Discovered in BradleY_et_al.i f2011) 

'Discovered in Trenti et al. {20lJ) 

sin order listed in table, discovered as P34.Z.6106, P34.Z.4809 in lWilkins et al.l JMTiT) 

•^In order listed in table, identified as ERS.z.45856 and ERS.z. 87209 in I Wilkins eTall l(20Tll 1 



broader than the rece nt vifork of lOno et al.l ()2011|) and 
iPentericci et al.l ()20lH ). In similar spectroscopic follow- 
up campaigns, they target brighter dropouts, primarily 
with Muv < -21.4, and -21.75 <Muv < -20.0, re- 
spectively. 

Remarkably, from the new Keck exposures, we find 
very few convincing detections of line emission. Figure 2 
(bottom panel) shows a 2-D spectroscopic montage of 4 
sources for which line emission may be present, of which 
one case (HUDF09.1596) is marginal {2a) and the other 
(ERS 8290) lies outside the expected redshift range if it 
is Lya. 

The emission line apparently seen in ERS 8290 (a z'- 
band dropout) is detected at >5ct with a flux of 6.7 ± 0.8 
X 10^^^ erg cm~2 s"^ at A7644 ± 2 A. It also exhibits 
the asymmetric profile characteristic of Lya in the ID 
extraction, but this would place the object at z = 5.286, 
quite discrepant from our photometric estimate of z = 
6.52. However, upon examining the positioning of our 
LRIS slits more carefully (bottom panel of Figure 2), we 
find there is a faint I^-drop candidate with i' = 27.5 and 
z' — 27.9 only 0.4" away which would have been at least 
partially visible through our slit during the exposure. 
After subtracting the line flux from the i'-band photom- 
etry, we determine a photometric redshift of z = 4.91, in 
reasonable agreement with the Keck spectroscopy, par- 
ticularly given a greater line flux (as is likely given the 
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Fig. 



1. — Distribution of rest-frame UV absolute magnitudes, 
MjjY, for i'-drop sources discussed in Paper II with z ~ 5.5-6.3 (top 
panel) compared with those for the present survey of i' and z'-drops 
at z > 6.3 (lower panel). There is an additional dropout, A2261.1, 
not shown on this histogram at Myy ~ —16. Dark shading in 
the lower panel refers to sources selected on the basis of WFC3/IR 
imaging in the Keck campaign, (Table 1); light shading r efers to 
additi onal data drawn from the VLT campaign of .Fontana et al.l 
120TO[) . 
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object's poor positioning in the slit) would increase the 
photometric redshift estimate. The resolution of the con- 
fusion arising from these two proximate sources empha- 
sizes again the prominence of Lya emission in low lumi- 
nosity 2: ~ 5 — 6 sources (Papers I and II). 

The two satisfactory detections refer to emission at 
A9045 ±2 A for ERS 8496 in the LRIS mask and emis- 
sion at A9780 ± 4 A in the NIRSPEC exposure of the 
lensed source Abell 1703 jD6. Both objects are detected 
at > 5(7 in our final exposures. In our ID extraction 
of ERS 8496, the emission line has a FWHM of 9 ± 1 
A, and displays an asymmetric profile with a steeply ris- 
ing blue edge and slowly decaying red tail, characteristic 
features of Lya at high redshift. Because our spectral 
resolution is significantly lower (6.5 A for NIRSPEC ver- 
sus 4.1 A for LRIS) in our spectrum of A1703 jD6, so we 
are unable to determine any line profile information. The 
emission feature is seen independently in coadditions on 
two separate nights, indicating its reality. 

If both lines are Lya, the implied spectroscopic red- 
shifts for ERS 8496 and Abell 1703^06 are z = 6.441 ± 
0.002 and 7.045 ± 0.004, respectively, in excellent agree- 
ment with our photometric predictio ns of 6.52 and that of 
7.0 derived bv iBradlev et al.l ()20lH ). The measured line 
fluxes for the two objects are 9.1 ± 1.4 and 28.4 ± 5.3 
X 10^^^ erg cm~^ s~ . We then assume a spectral slope 
of /3 = -2, which is chara cteristic of galaxies at this red- 
shift (jDunlop et al.ll201ll ) , though there n iay be evidence 
for ste eper slopes at sub-L* luminosities (jBouwens et al.l 
l2010al) . Under this assumption, taking the magnitude 
from the first filter in which the object is detected (1o98A/ 
for ERS 8496, and J125 for 1703 jD6), we calculate EWs 
of 69 ± 10 and 65 ± 12 A, respectively. Because our ob- 
jects have additional coverage longward of the detection 
filter, we can also compute a value for /3, and extrapolate 
to find the continuum fiux at \r est = 1216 A. Using the 
formulae of iDunlop eFall ()201lD . we find /3 = -2.39 ± 
0.55, and -2.44 ± 0.64. When computing EWs using this 
method, we obtain 67 ± 11 and 59 ± 12 A, respectively. 

3. ANALYSIS 

Our approach in this paper is to compare the rate of 
occurrence of Lya in our new 6.3 < z < 8.8 sample with 
that expected from our reference sample of i'-drops with 
5.5 < z < 6.3 drawn from Paper II (Figure 1). In both 
Papers I and II we showed that the rest-frame EW dis- 
tribution is a function of rest-frame UV absolute magni- 
tude, Muvj E^nd thus we additionally take this luminosity 
dependence into account. We estimate the luminosities 
of all our sources in Table 1 from their photometric red- 
shift and incorporat e the lensing m agnification fi for our 
lensed sources from IBradlev et aT (2011) for Abell 1703 
and from the mass model of Richard ct al. (2010) for 
MS0451-03 and Abell 2261. For our baseline Lya EW 
distribution, we use the data from of Paper II at z ~ 6, 
separated into high and low luminosity regimes. 

The fraction of emitters within each bin of EW > 25 
A is taken directly from Fig. 2 of Paper II. We set the 
slope of the distribution within an EW bin equal to the 
slope between the two lowest bins in Paper II, 25 A < 
EW < 55 A, and 55 A < EW < 85 A. This slope is 
equal to dp(EW)/dEW = —0.0030 for the lower luminos- 
ity sample (-20.25 <Muv < -18.75), and -0.0017 for 
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Fig. 2. — Montage of Lyo emission detected from 4 sources in 
our Keck survey, along with boxcar-extracted ID spectrum. Wave- 
length ranges contaminated by strong skylines are shaded in grey 
in the ID extraction. The top row shows the two robust detections 
of ERS 8496 and A1703_zD6 at z = 6.441 and z=7.045 respec- 
tively. The bottom row shows a marginal detection at z=6.905 for 
HUDF09_1596 and a likely Lya line at z = 5.286 arising from a 
serendipitous V-drop close to ERS 8290 as illustrated in the bot- 
tom left slit image. The bottom right panel shows the photometric 
redshift distribution for this serendipitous V-drop, with a verti- 
cal line indicating the observed spectroscopic redshift (see text for 
further details). 



the higher luminosity sample (—21.75 kMjjv < —20.25). 
To create the probability distribution for galaxies with 
EW less than 25 A, we extrapolate to EW = A using 
this slope, and assign the remaining fraction of galaxies 
as non-emitters. In Papers I and II we also showed the 
fraction of emitters is a function of redshift, rising sig- 
nificantly for lower luminosity sources over 4 < z < 6, 
most likely as a result of reduced dust extinction in the 
early Lyman break population. Therefore, as discussed 
in Paper II, we have also used a projected rest-frame EW 
distribution at z ~ 7 assuming this evolutionary trend 
continues beyond z ~ 6. 

Two key factors enter into the calculation of the visibil- 
ity of line emission in a ground-based survey. Firstly, for 
any target with a particular photometric redshift like- 
lihood function p{z), it may be that the spectral re- 
gion surveyed by LRIS or NIRSPEC does not completely 
cover the expected wavelength range where Lya might be 
present. Secondly, the EW limit for Lya emission will 
be a highly non-uniform function of wavelength due to 
the mitigating effect of night sky emission. Provided the 
photometric redshift solution we derive is robust, we can 
estimate both factors and hence derive the likelihood of 
seeing Lya for each of our 26 sources assuming the rele- 
vant wavelength range studied and the exposure time se- 
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cured, if the particular source of a given Muy has a EW 
distribution drawn from the sample with 5.5 < z < 6.3. 

In the case of those sources for which the wavelength 
range searched does not fully sample the extended p{z), 
we reduce the detection likelihood by the fraction of 
the integrated p{z) that lies outside our search range. 
For each target, we determine its redshift probability 
function p(z) using th e photometric redshift code EAZY 
(jBrammer et al.|[2008() . To determine the varying visibil- 
ity function within our search range, we first estimate the 
the noise within an aperture encompassing the expected 
profile of the line, assuming an emission line width of 10 
A FWHM which is typical of those detected in Paper 
II. Figure 3 illustrates the 5a EW limit as a function of 
wavelength (and Lya redshift) for most sources in our 
12 hour LRIS exposure (6.3 < z < 7.2) and a typical 
source studied with NIRSPEC (6.8 < z < 8.2) during a 
5 hour exposure. We note that although our survey spans 
a large range in redshift (6.3 < z < 8.8), the p{z) dis- 
tributions for individual galaxies typically span a much 
smaller range. The average la redshift confidence inter- 
val for sources in our sample is only Az = 0.43. 

Since the NIRSPEC exposures were usually single- 
object exposures, the limits vary from source to source 
depending on the conditions and exposure times. We 
then apply a completeness correction to account for the 
fact that an emission line may fall in a noise trough and 
lie undetected, despite having an intrinsic fiux above the 
5ct limit. To estimate this completeness correction, we 
follow the methods discussed in detail in Paper II, where 
we simulate the addition and recovery of fake line emis- 
sion in our actual spectra, again assuming a FWHM of 
10 A. As the absolute limits vary from source to source 
our Keck survey is not complete to a fixed EW limit but, 
provided the limits are well-understood for each source, 
we can readily estimate the proba bility of seeing Lyct in 
our exposures. In the case of the iFontana et al.l ()2010l ) 
F0RS2 survey we estimated the night sky emission from 
our own LRIS exposures normalizing the limits from nu- 
merical data supplied in that paper. 

The above simulations can be used to verify that our 
Keck survey is well-placed to search for Lya emission. 
Out o f the combined 26 t argets from our survey and 
that of lFontana et al.l (|2010D , 24 are covered spectroscop- 
ically over more than half the integrated probability of 
their photometric redshift distribution, and 17 are cov- 
ered over 95% of the range. Additionally, we are able 
to determine the fraction of our spectra occulted by OH 
sky emission. For example, for a J=27 galaxy in one of 
our LRIS exposures, we are sensitive to lines with EW > 
30 A over 70% of our usable spectral range (see Figure 
3). Similarly, for a 5 hr NIRSPEC exposure of a J= 27 
galaxy, we are sensitive to lines with EW > 55 A over 
49% our spectral range. 

The results of the simulations are shown in Figure 4 
where, depending on whether we adopt the EW distri- 
bution observed at 5.5 < z < 6.3 or that extrapolated 
to 2; ~ 7 in Paper II assuming continuity in redshift- 
dependent increase in line emission seen over 4 < z < 6, 
we would expect to recover 7-8 emission lines in the com- 
bined Keck and VLT surveys. In contrast, we have only 
2 robust detections (both in the Keck sample, Figure 2) 
and at most 4 including the marginal candidate discussed 




Fig. 3. — Sensitivity limits to Lyo emission in our new Keck 
spectroscopic campaign: The panels show the 5cr limiting EW cal- 
culated for a typical source studied in our 12 hour multi-slit LRIS 
exposure (top) and an example 5 hour long slit NIRSPEC observa- 
tion (bottom). The limits vary from source to source depending on 
the continuum brightness and the exposure times. An additional 
completeness correction is taken into account by adding and at- 
tempting to recover fake emission lines with fluxes equal to the Scr 
flux limit at the wavelength of insertion. See text and Paper II for 
more detailed discussion. 



bv IFontana et al.l (j2010( ) and HUDF09_1596 at z = 6.905 
shown in Figure 2. Assuming all the targeted sources 
are at z > 6.3, given our previously mentioned assump- 
tions, our results reject the input EW distributions at the 
99.3% level of significance (91.4% if the two marginal de- 
tections are included). 

We can display the significance of this downturn with 
increasing redshift in the terms of the fraction of Lya 
emission seen in Lyman break galaxies, X(Lya) as in 
Paper I. The difficulty we face in creating such a fig- 
ure is the non-uniform EW limit across the various tar- 
gets in the Keck and VLT campaigns, in contrast to the 
more straightforward uniform search we undertook with 
DEIMOS at 4 < z < 6. To account for this, we assume a 
simple model in which Lya emission is transmitted with- 
out IGM absorption for a fraction / of galaxies, while it 
is fully extinguished by the IGM for a fraction of galaxies 
(1-/). We assume that / = 1 at redshifts belo w 6, where 
the u niverse is believed to be highly ionized (jFan et al.l 
I200(]f ). and that / is independent of the intrinsic EW of 
a Lya emission line. We caution that an interpretation 
in terms of absolute values of / is premature, as there is 
still some debat e on w hether the IGM is fully ionized at 
z ~ 6 Jvlcsingei^ (j2010[ ) , but emphasize that our value of 
/ at z ~ 7 is computed relative to the value assumed at 
z ~ 6. Additionally, with the increased fraction of emit- 
ters in our z '^ 6 sample from Paper II, we do not see 
any evidence for a decrease in / prior to z ~^ 6, though 
we cannot rule it out. 

It is important to note that our / is differs from /^f^^", 
commonly defined in the liter ature as the total escape 
fraction of Lya photons (e.g., iHaves et al.ll201ll) . f^J^^" 
represents the total transmission of Lya, accounting for 
both attenuation of photons within the galaxy by mech- 
anisms such as dust, as well as any attenuation by the 
IGM. Our definition of / is only intended to account for 
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Fig. 4. — The expected number of detected Lya emission lines 
with greater than or equal to Scr significance in the combined Keck 
and VLT survey of 26 sources. The blue histogram shows the like- 
lihood function for 10,000 Monte Carlo realizations assuming the 
intrinsic line emission properties follow the luminosity dependence 
seen in our 5.5 < z < 6.3 i'-drop sample (Figure 1 (top). The open 
histogram shows the expectation if the fraction of line emitters con- 
tinues to increase with redshift at the rate described in Paper II. 
Vertical lines show the recovered number of emitters (robust and 
maximal including marginal detections in both the Keck and VLT 
surveys) . 



any downturn in the fraction of LBGs which show observ- 
able Lya emission from the z = 6 (or z = 1) extrapolated 
EW distributions from Paper II, and represents an IGM 
extinction averaged over the entire population. 

To compute the most likely value of /, we undertake 
Monte Carlo simulations using the previously described 
EW distributions, but with / now added as a free param- 
eter. We vary / from to 1 in steps of 0.01, and compute 
N=1000 simulations for each step. We can then calcu- 
late the probability distribution for / given our No(,s=2 
confirmed sources using Bayes' theorem: 



p{f\Nobs 



pjNobs = 2|/)p(/) 



(1) 



Here, p{f) is the prior probability for /, which we take 
to be uniform for < / < 1, and p{Nobs = 2|/) is the 
probability, drawn from our Monte Carlo simulations, 
that we would find ^obs = 2 sources for a given value 
of /. Assuming that the intrinsic EW distribution for 
our observed sources is that of Paper II at z = 6, we 
find / = 0.45 ± 0.20, while using the z = 7 extrapolated 
distribution yields / — 0.34^Q^g. In the Figure 5, we 
plot the value of X(Lya) in the same luminosity bins of 
Paper II, as predicted by our best fit values of /. 

We stress that this figure is intended to serve as a con- 
tinuation of the visualization provided in Papers I and II, 
rather than a statistiscal result of our study. Due to our 
strongly varying limiting EW sensitivity (as a function of 
both wavelength and object magnitude), choosing a fixed 
EW limit will exclude a non-negligible fraction of useful 
data from our analysis. Our Monte Carlo simulations are 
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Fig. 5. — The redshift-dependent fraction of color-selected Ly- 
man break galaxies that reveal Lya in emission, X(Lya), adjusted 
as discussed in the text to approximate one within a similar lumi- 
nosity range with a rest-frame EW in excess of 25 A. Data points 
for tlie galaxies with —21.75 < Muv < —20.25 are displaced by 
-1-0.1 in redshift for clarity. Data over 4 < 2; < 6 is from Paper I 
and Paper II, and new estimates beyond z > 6.3 ar e derived from 
the p resent paper, including sources discussed by IFontana et al.l 
II2OIOI ). The curves shown represent the aggregate redshift proba- 
bility distributions for our sources in the z ~ 6 bin (black), and the 
z ~ 7 bin (blue); probability distributions for individual sources 
are typically much sharper. 



able to utilize the full data set, simulating whether we 
would have likely seen a line even when our EW limits 
are above the fixed thresholds used in Fig. 5, and thus 
represent the major stat istical result of this stu dy. 

Using the models of iMcQuinn et al.l ()2007() to pre- 
dict what global neutral hydrogen fraction, Xhi would 
be required to account for this decline, we find Xhi — 
0.44, and Xh t — '^^' respectively. The models of 
iDiikstra et al.l (|201lD . which provide a more comprehen- 
sive treatment of Lya radiative transfer through out- 
flows, result in an increased value for Xjji in both cases. 



4. DISCUSSION 

Although we consider the most likely explanation for 
our observed decrease in the number of LBGs which show 
observable Lya emission to be an increase with redshift 
in the neutral fraction of the IGM, it is important to re- 
member our assumptions. Foremost we have assumed 
that all of our 26 targets have true redshifts beyond 
z ~ 6.3. Should there be low redshift interlopers or 
Galactic stars in our new sample, we will overestimate 
the decline in the Lya fraction. Secondly, we have as- 
sumed the DEIMOS spectra from Paper II constitute a 
representative sample for calculating the expected EW 
distribution for 6.3 < z < 8.2. Although the uncertain- 
ties here are not as great, we plan further studies with 
DEIMOS to increase the statistical sample of 5.5 < z < 
6.3 LBGs. 

Of course our observed decrease in the Lya fraction 
could also be attributed to an increased opacity arising 
from dust within the LBGs. However, given the blue 
UV continuum slopes observed for galaxies with z > 6.3 
(jBouwens et al.., 2010a : .Dunlop et al.. i2011.), we consider 
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this explanation unlikely. 

Our diagnosis of a possible increase in the neutral hy- 
drogen fr action beyond z c± 6.3 is supported by the earlier 
study of iFontana et al.l ()2010f ). They found 1 marginal 
candidate out of 7 targets whereas we find 2 robust 
and 1 marginal cases out of our 19 targets spanning a 
larger luminosity and redshift range. Our conclusion is 
also supported by LAE s studies at z = 5 .7 an d 6.5 by 
lOuchi et all (|2010l ) and iKashikawa et al.l poTTI ). Com- 
pared to z = 5.7, their LAE sample at z = 6.5 displays 
systematically lower EWs for Lya. They also derive little 
evolution in the rest UV luminosity function for LAEs, 
but a decrease in the Lya LF, which could be explained 
by an increase in X/f/. O ur derived values of Xg j are 
slightly higher than that of IKashikawa et aTl (|201lD . per- 
haps consistent with our surve y probing to highe r red- 
shifts than their z = 6.5 LAEs. Hav es et all ([2011) have 
recently compiled results from numerous Lya and UV 
luminosity function studies to derive a volumetrically av- 
eraged Lya escape fraction, and find very similar results. 
Their derived Lya escape fraction steadily increases with 
redshift below = 6, then tentatively drops off at higher 
redshifts. 

Very recently, lOno et aTl (|2011[ ) report the convinc- 
ing detection of Lya emission in a small fraction (3/11) 
of LBGs that, by virtue of their selection using Subaru 
imaging, are more luminous (Myv < ~21) than most 
of the objects considered here. Such a complementary 
campaign targeting luminous LBGs selected from larger 
volumes will provide further insight into whether reion- 
ization is responsible for the declining fraction of line 
emission. 

We note that our measured decrease in the fraction of 
LBGs with strong Lya potentially agrees with the result 
of Cowic et al. (2011). Although they argue against any 
evidence for reionization at z = 6.5, they find that ~ 24% 
of galaxies at this redshift show strong Lya emission, 
comparable to the fraction we detect in this work, spread 
across a larger redshift range. 

With the new generation of multi-object, near infrared 
spectrographs, such as MOSFIRE, set to come online 
soon, the prospects for this field are bright. In addition 
to the significant multiplexing advantage, the increased 
sensitivity of these detectors will allow us to probe the 
lower luminosity ranges at z > 6.5 to EW limits compa- 
rable to those in Paper II between sky lines. Having such 
a statistical sample is key for allowing the quantification 
of any change in the hydrogen neutral fraction. 
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